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Prethermalization of the equation of state and the kinetic temperature to their equihbrium values 
occurs on time scales dramatically shorter than the thermal equilibration time. This is a crucial 
ingredient for the understanding of collisions of heavy nuclei or other nonequilibrium phenomena 
in complex quantum and classical many body systems. We also compare the chemical equilibration 
time with other characteristic time scales. 
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PACS numbers: 11.10.Wx,12.38.Mh,05.70.Ln 

Prethermalization is a universal far-from-equilibrium 
phenomenon which describes the very rapid establish- 
ment of an almost constant ratio of pressure over energy 
density (equation of state) , as well as a kinetic tempera- 
ture based on average kinetic energy. The phenomenon 
occurs on time scales dramatically shorter than the ther- 
mal equilibration time. As a consequence, prethermal- 
ized quantities approximately take on their final thermal 
values already at a time when the occupation numbers 
of individual momentum modes still show strong devi- 
ations from the late-time Bose-Einstein or Fermi-Dirac 
distribution. 

The abundance of experimental data on matter in ex- 
treme conditions from relativistic heavy-ion collision ex- 
periments, as well as applications in astrophysics and cos- 
mology urge a quantitative understanding of nonequilib- 
rium dynamics in quantum field theories. Collision ex- 
periments seem to indicate early thermalization whereas 
the present theoretical understanding of QCD suggests 
a much longer thermal equilibration time. For example, 
the successful application of hydrodynamics already less 
than 1 fm after the collision is so far unexplained from 
theory 

To resolve these questions, it is important to under- 
stand to what "degree" thermalization is required to ex- 
plain the observations. Different quantities effectively 
thermalize on different time scales and a complete ther- 
malization of all quantities may not be necessary. For in- 
stance, an approximately time-independent equation of 
state p = p(e), characterized by an almost fixed rela- 
tion between pressure p and energy density e, may form 
very early — even though the system is still far from 
equilibrium! Such an almost constant equation of state 
is a crucial ingredient for the use of efficient hydrody- 
namic descriptions, since it is needed to close the system 
of equations obtained from the conservation of the energy 
momentum tensor. 

The initial stages of a collision require to consider 
quantum fields in extreme nonequilibrium situations. 
Connecting this far-from-cquilibrium dynamics at early 
times with the approach to thermal equilibrium at late 
times is a challenge for theory. Achieving this goal is cru- 



cial for a comparison between the time scales of prether- 
malization and thermal equilibration, and we investigate 
this question here quantitatively in an effective quark- 
meson model. Approaches based on small deviations 
from equilibrium, or on a sufficient homogeneity in time 
underlying kinetic descriptions, are not applicable to de- 
scribe the "link" between the initial and the late-time 
behavior. Classical field theory approximations are ex- 
pected to be valid for not too late times, but cannot de- 
termine the relevant time scale for the approach to quan- 
tum thermal equilibrium. Recently, it has been demon- 
strated IJi 1^ 1^ that far-from-equilibrium dynamics 
as well as subsequent thermalization of quantum fields 
can be described using efficient functional integral tech- 
niques. The description includes direct scattering as well 
as off-shell and memory effects. This is crucial to es- 
tablish the different time scales for a loss of memory of 
initial conditions for certain "bulk quantities" as com- 
pared to "mode quantities" characterizing the evolution 
of individual momentum modes. 

The observation that the nonequilibrium system looses 
a major part of the memory of the detailed initial con- 
ditions on a very short time scale is a robust feature of 
classical as well as quantum field theories. It has been 
observed that approximate "equipartition" between 
bulk kinetic and potential energy occurs very rapidly. 
Below, this will be the basis of our definition of a "ki- 
netic temperature" Tkin- We will see that the equation of 
state becomes almost constant at the same "prethermal- 
ization time" ipt. A rapid approach to a slow evolution 
of the equation of state in classical field theories has also 
been observed with expanding space-time geometries |^ . 
The fast loss of memory for these quantities is based on 
the phenomenon of "dephasing" , which is independent 
of the interaction details. In contrast, "mode tempera- 
tures" (to be defined below) for individual momentum 
modes loose only part of the initial condition details on 
a somewhat longer time scale tdamp which depends on 
the interaction strength llSlj . Still, tdamp is much smaller 
than the true thermal equilibration time toq which de- 
scribes the universal rate of approach to the equilibrium 
values for all relevant correlation functions y, lol - An 
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FIG. 1: Fermion occupation number ■nS^^{t\p) for three dif- 
ferent momentum modes as a function of time. The evolution 
is shown for two different initial conditions with same energy 
density. The long-time behavior is shown on a logarithmic 
scale for t > 30 m~^. 



even substantially larger separation of scales is observed 
in classical field theories 0, 0, 0| as compared Q to 
the corresponding quantum theories. 

In this letter we consider the nonequilibrium evolution 
of quantum fields for a low-energy quark-meson model, 
which takes into account two quark flavors with a Yukawa 
coupling ~ /i to a scalar cr-field and a triplet of pseu- 
doscalar pions, vr. The theory corresponds to the well- 
known "linear cr-model" , which incorporates the chiral 
symmetries of massless two-flavor QCD. The classical ac- 
tion reads 

S ^ j d'*a;|V'i^V' + \ [^^c^S'^ct + a^TrS'^Tf] 

+ [fj -f Z75f7r] -0 - + 7r^)| . (1) 

We consider a quartic scalar self-interaction V{(j'^+tt'^) — 
ml + 7r2) /2 + A (a^ + ^2)^ /{m]) with Nf = 2. 
The employed couplings are taken to be of order one, and 
if not stated otherwise h = X = 1. We emphasize that the 
main results of this letter about prethermalization are in- 
dependent of the detailed values of the couplings. Here 
we use the two-particle irreducible (2PI) effective action 
to two- loop order . In Ref . it has been shown that 
this approximation can be used to study the far-from- 
equilibrium dynamics as well as the late-time approach 
to quantum thermal equilibrium. The dynamics is solved 
numerically without further approximations (cf. Ref. 
for calculational details). All quantities will be given in 
units of the scalar thermal mass m j^^l • 

Thermalization: Nonequilibrium dynamics requires 
the specification of an initial state. A crucial question of 
thermalization is how quickly the nonequilibrium system 
effectively looses the details about the initial conditions, 
and what are the characteristic stages of a partial loss 
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FIG. 2; Fermion and scalar mode temperatures Ty'''\t) as a 
function of momentum p for various times. 



of information. Thermal equilibrium keeps no memory 
about the time history except for the values of a few con- 
served charges. In Fig.Qlwe show the effective occupation 
number density of fermion momentum modes, vS^Ht\ p], 
as a function of time for three different momenta |2l|. 
The plot shows two runs denoted as (A) and (B) with 
different initial conditions but same energy density. Run 
(A) exhibits a high initial particle number density in a 
narrow momentum range around ±p. This situation is 
reminiscent of two colliding wave packets with opposite 
and equal momentum. We emphasize, however, that we 
are considering a spatially homogeneous and isotropic en- 
semble with a vanishing net charge density. For run (B) 
an initial particle number density is employed which is 
closer to a thermal distribution. 

One observes that for a given momentum the mode 
numbers of run (A) and (B) approach each other at early 
times. The characteristic time scale for this approach 
is well described by the damping time idamp(p) [13 ■ 
Irrespective of the initial distributions (A) or (B), we 



find (for p/m ~ 1) t 
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is) 



28 m- 



damp 

for scalars. 



25 m ^ for fermions and 
In contrast to the initial 



rapid changes, one observes a rather slow or "quasista- 
tionary" subsequent evolution. The equilibration time 
icq — 95 m~^ is substantially larger than tdamp and is ap- 
proximately the same for fermions and scalars. Thermal 
equilibration is a collective phenomenon which is, in par- 
ticular, rather independent of the momentum. In sum- 
mary, mode quantities such as effective particle number 
distribution functions show a characteristic two-stage loss 
of initial conditions: after the damping time scale much 
of the details about the initial conditions are effectively 
lost. However, the system is still far from equilibrium 
and thermalization happens on a much larger time scale. 

We define mode temperatures Tp^'^^ (t) by equating the 
mode particle numbers ri'p'^\t) with a time and momen- 
tum dependent Bose-Einstein or Fermi-Dirac distribu- 
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FIG. 3: The ratio of pressure over energy density ui as a func- 
tion of time. The inset shows the early stages for two different 
couplings and demonstrates that the prethermalization time 
is independent of the interaction details. 



tion, respectively 



np{t) = [exp (ujp{t)/Tp{t))±lY 



(2) 



This definition is a quantum mechanical version of its 
classical counterpart as defined by the squared "gen- 
eralized velocities" .7]. In thermal equilibrium with 
ujp ~ ^/^PTJP and Tp = T^q Eq. @ yields the fa- 
miliar occupation numbers (/k = 0). Here the mode fre- 
quency ajp^'^\t) is determined by the peak of the spec- 
tral function for given time and momentum, as detailed 
in Ref. . In Fig. |21 we show the fermion and scalar 
mode temperature as a function of momentum for var- 
ious times t 3> idamp- One observes that at late times, 
when thermal equilibrium is approached, all fermion and 
scalar mode temperatures become constant and agree: 



'{t) = Tp^\t) — Tcq. In contrast, there are sizeable 
deviations from the thermal result even for times consid- 
erably larger than the characteristic damping time. 

Kinetic prethermalization: In contrast to the rather 
long thermalization time, prethermalization sets in ex- 
tremely rapidly. In Fig. |21 we show the ratio of pres- 
sure over energy density, w = p/e, as a. function of time. 
One observes that an almost time-independent equation 
of state builds up very early, even though the system is 
still far from equilibrium! The prethermalization time 
tpt is here of the order of the characteristic inverse mass 



scale 



This is a typical consequence of the loss of 



phase information by summing over oscillating functions 
with a sufficiently dense frequency spectrum. In order 
to see that this phenomenon is not related to scattering 
or to the strength of the interaction, we compare with a 
smaller coupling in the inset and observe good agreement 
of both curves. The dephasing phenomenon is unrelated 
to the scattering-driven process of thermalization. 

Given an equation of state, the question arises whether 
there exists a suitable definition of a global kinetic tem- 
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FIG. 4: Chemical temperatures for scalars (upper curves) and 
fermions (lower curves) for different values of the coupling h 
and Tcq. We also show the kinetic temperature rkin(t) (solid 
line) , which prethermalizes on a very short time scale as com- 
pared to chemical equilibration. 



perature Tkin. In contrast to a mode quantity such as 
Tp (t) , a temperature measure which averages over all mo- 
mentum modes may prethermalize. Building on the clas- 
sical association of temperature with the mean kinetic 
energy per degree of freedom, we use here a definition 
based on the total kinetic energy £'kin(i): 



rkin(t) =-Bkin(t)/c, 



cq ■ 



(3) 



Here the extensive dimensionlcss proportionality con- 
stant Ccq = -Ekin rn/Tpn is giveu solely in terms of equilib- 
rium quantities [23 • Since total energy is conserved, the 
time scale when "equipartition" is reached (i.e. E'kin/^^ 
is approximately constant) also corresponds to a time- 
independent kinetic temperature. The latter equals the 
equilibrium temperature Teq if £'kin/£' has reached the 
thermal value. 

The solid line of Fig.Qlshows Tkin(i) normalized to the 
equilibrium temperature (for Tcq/m — 1). One observes 
that an almost time-independent kinetic temperature is 
established after the short-time scale tpt ~ m~^. The 
time evolution of bulk quantities such as the ratio of 
pressure over energy density w, or the kinetic temper- 
ature Tkin, are dominated by a single short-time scale. 
These quantities approximately converge to the thermal 
equilibrium values already at early times and can be used 
for an efficient "quasi-thermal" description in a far-from- 
equilibrium situation! 

Chemical equilibration: In thermal equilibrium the rel- 
ative particle numbers of different species are fixed in 
terms of temperature and particle masses. A system 
has chemically equilibrated if these ratios are reached, 
as observed for the hadron yields in heavy ion colli- 
sions 01 . Obviously, the chemical equilibration time 
tch will depend on details of the particle number chang- 
ing interactions in a given model and tch £ icq- In 
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our model we can study the ratio between the num- 
bers of fermions and scalars. For this purpose we in- 
troduce the chemical temperatures Tl:^'^\t) by equating 
the integrated number density of each species, n^-'"''*^ (t) — 
g{f,s) ^ A^p/{2TTfni!'''\t), with the integrated Bose- 
Einstein/Fermi-Dirac form of distributions: 

/•OO 

n{t) = / dpp" [exp (a;p(i)/r,h(t)) ± 1]"' . (4) 

Here g^^'^ — 8 counts the number of fermions and g'*^ = 4 
for the scalars. 

The time evolution of the ratios T^^^\t) /T^t^ is shown 
in Fig. ^ for different values of the coupling constant h 
and the equilibrium temperature T^q. One observes that 

chemical equilibration with T^^\t) — T^l\t) does not 
happen on the prethermalization time scale, in contrast 
to the behavior of T]iin{t). Being bulk quantities, the 
scalar and fermion chemical temperatures can approach 
each other rather quickly at first. Subsequently, a slow 
evolution towards equilibrium sets in. For the late-time 
chemical equilibration we find for our model tch — icq- 
However, the deviation from the thermal result can be- 
come relatively small already for times t <C tcq. 

Let us finally consider our findings in view of colli- 
sions of heavy nuclei and try to estimate the prethermal- 
ization time. Actually, tpt is rather independent of the 
details of the model like particle content, values of cou- 
plings etc. It mainly reflects a characteristic frequency 
of the initial oscillations. If the "temperature" (i.e. av- 
erage kinetic energy per mode) sets the relevant scale 
one expects Ttpt = const. (For low T the scale will 
be replaced by the mass.) For our model we indeed find 
Ttpt —2 — 2.5 We expect such a relation with a sim- 
ilar constant to hold for the quark-gluon state very soon 
after the coUision For T > 400 - 500 MeV we obtain 
a very short prethermalization time tpt of somewhat less 
than 1 fm. This is consistent with very early hydrody- 
namic behavior [2^. In QCD the equilibrium equation 
of state shows no strong temperature dependence above 
the critical temperature [ig , and can therefore adapt 
easily as the temperature decreases. After the transition 
w will only re-adjust somewhat to the equilibrium value 
relevant for a hot hadron gas, typically on a time scale of 
a few fm (for T ~ 175 MeV). The chemical equilibration 
time tch depends on the production rate for multi-strange 
hadrons [l3- From tpt <C tch and Fig. ^ we conclude that 
once the chemical temperatures for the various differ- 
ent species are equal the relevant chemical temperature 
Tch coincides with Tkin and defines a universal tempera- 
ture. Comparison with the critical temperature in equi- 
librium is therefore meaningful - an approximate equality 
Tch — Tc has been advocated ^3 ~ such that chemical 
freeze out can, in principle, be used to measure Tc. 
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